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CAPTIONS 

1.  The  SSJ/3  instrument  and  associated  electronic 
boards.  The  rectangular  slits  are  entrance  apertures 
for  electrons.  Behind  the  small  circular  hole  is 

a UV  detector  that  causes  the  channeltron  acceleration 
voltage  to  be  turned  off  when  viewing  the  sun. 

2.  Detailed  phtotgraph  of  the  SSJ/3  head  assembly  board, 
a h"  block  is  also  shown  for  comparison.  The 
narrowly  separated  plates  select  the  high-energy 
electrons.  The  widely  separated  plates  select  the 
lower  energy  electrons. 

3.  Detailed  schematic  sketch  of  a detector  assembly. 

Also  given  are  the  relevant  dimensions  for  each 
set  of  plates  . 

4.  Schematic  of  the  electronic  logic  used  in  SSJ/3. 

5.  Comparison  between  measured  and  Monte-Carlo  results 
for  a normally  incident  beam.  Energies  are  off-set 
for  ease  of  display. 

6.  Angular  response  of  the  high-energy  plates  in  the 
R-plane  (a) . 

7.  Angular  response  of  the  low  energy  plates  in  the 
R-plane  (a). 

8.  Angular  response  of  the  high-energy  plates  in  the 
plane  parallel  to  the  plates  fS). 

9.  Angular  response  of  the  low  energy  plates  in  the 
plane  parallel  to  the  plates  ('3). 

10.  Theoretical  response  of  the  high-energy  plates  to 
an  isotropic  incident  flu.x.  Use  Table  2 for 
channel  8 mid-point  energy.  Final  normalization 
is  also  shown. 


11.  Theoretical  response  of  the  low  energ\’  plates  to 
an  if  pic  incident  flux.  Use  Table  2 for 
channel  9 mid-point  energy.  Final  normalization 
is  also  shown. 
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I.  INTRODUCTION 


The  calibration  of  the  SSJ/3  sensor  was  carried  out 
using  both  experimental  and  theoretical  techniques.  Experi- 
mental measurements  were  compared  with  predictions  from  a 
Monte  Carlo  particle  tracing  code  and  from  analytic  approxi- 
mations . 

In  Section  II  instrument  design  is  discussed.  Section  III 
described  techniques  used  in  calibration  while  Section  IV 
describes  the  logic  used  to  obtain  the  final  normalization 
values . 


5 


..  . . 


II  . INSTRUMENT  DESIGN 

Introduction . The  SSJ/3  instrument  package  is  shown  in 
Figures  1 and  2.  Figure  1 shows  the  instrument  case  and  the 
associated  electronic  boards.  The  apertures  on  the  case  are 
directly  in  front  of  the  particle  collimators  that  are  shown 
in  the  bottom  of  Figure  2.  Behind  the  collimators  are  curved 
cylindrical  plates.  The  plate  separation  and  the  applied 
voltage  determine  the  energy  band  of  the  transmitted  electrons 
The  electrons,  after  passing  through  the  exit  collimator,  are 
detected  by  the  C-shaped  Channeltron  Electron  Multiplier 
detectors  that  are  also  shown  in  Figure  2. 

Head  Assembly.  Electrostatic  analyzers  (ESA's)  select 
charged  particles  by  applying  a voltage  between  two  con- 
centric plates.  The  SSJ/3  instrument  uses  two  sets  of 
cylindrical  concentric  plates  as  shown  in  Figure  2.  Charged 
particles  of  the  proper  sign  and  energy  have  trajectories 
that  are  almost  parallel  to  the  plate  surfaces.  Particles  of 
greater  or  lesser  energy  impinge  on  the  plate  surfaces.  The 
mean  energy  E can  be  easily  derived  by  requiring  a balance 
between  electrostatic  and  centrifugal  accelerations.  For 
plates  of  radii  rj  and  T2  (^2^  tj)  and  an  applied  voltage 
difference  of  V,  E is  given  by 


E = 


eV 


2 ln(r2/ri) 


(1) 


If  V is  fixed  then  the  detected  energy  is  determined  by  the 
ratio  of  r2  to  r^.  This  explains  the  larger  plate  separation, 
AR,  for  the  lower  energy  detector  assembly  as  tabulated  in 
Figure  3.  The  sector  angle  of  the  plates  and  their  separations, 
together  with  aperture  and  exit  slits,  determine  the 
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transmission  efficiency.  Additional  baffles  are  placed  in  the 
aperture  assembly  to  eliminate  sunlight  (UV)  scatter  into  the 
detectors . 

Behind  the  exit  aperture,  the  two  channeltrons  are 
mounted  as  shown  in  the  lower  righthand  corner  of  Figure  2. 
Electrons  impinging  on  the  cone-shaped  area  produce 
secondaries  which  cascade  through  the  channeltron  producing  a 
detectable  pulse  at  the  end.  This  pulse  is  then  processed  by 
the  instrument's  internal  electronics.  Channeltron  operating 
lifetime  is  extended  by  shutting  off  the  accelerating  potential 
whenever  the  instrument  is  looking  at  the  sun.  This  is  done 
by  a small  phototransistor  mounted  adjacent  to  the  instrument 
apertures.  See  Figure  1.  In  Figure  3 is  a schematic 
representation  of  the  SSJ/3  head  assemblies  and  gives  the 
relevant  dimensional  values. 

Electronic  Logic.  Figure  4 shows  the  electronic  logic. 

The  channeltron  outputs  (DET  L and  DET  H)  are  amplified  and 
then  counted  in  nine-bit  logarithmic  counters.  A complete 
16-channel  energy  spectrum  is  read  out  every  second.  Appendix 
A describes  in  detail  the  conversion  of  the  logarithmic 
counter  to  decimal  form.  The  same  programmable  power  supply 
is  used  for  biasing  both  sets  of  plates.  The  voltage  dif- 
ference for  each  set  of  plates  is  obtained  by  applying  t\' / 7 
to  each  plate. 


III . TECHNIQUES  USED  IN  CALIBRATION 

Experimental . Experimental  calibrations  were  made  using 
accelerated  electrons  from  a Tritium  source.  After  acceleration 
by  an  applied  voltage,  the  beam  was  collimated  so  that  a mono- 
energetic,  unidirectional  source  was  incident  on  the  instrument 
apertures.  A turntable  enclosed  in  a vacuum  system  provided 
mobility  for  angular  response  measurements.  The  beam  was 
found  to  have  an  energy  dispersion  which  was  negligible  in 
comparison  with  the  instrument's  response.  By  changing  the 
applied  voltage  various  energies  could  be  selected. 

Theoretical  Approximations.  Geometric  considerations 
lead  to  the  following  results:  1)  If  a and  b are  the  heights 
of  the  entrance  and  exit  apertures  of  the  collimator  and  is 
the  separation  then  between  them  such  that  ^ >>  a + b,  then 
the  angular  resolution,  AB,  is  to  a good  approximation  given 
by  ASfti  (a  + b)/£  . This  approximation  works  well  for  the 

direction  parallel  to  the  plates,  but  not  for  the  orthogonal 
component.  2)  The  geometric  factor  (cm^-sec)  for  a long 
rectangular  collimator  is  approximated  by 
G = ^.L  X A2 

I 

- area  of  collimator  entrance  aperture 
A2  = area  of  collimator  exit  aperture 
I = distance  between  apertures 

This  formula  provides  a rapid  means  of  obtaining  a "ballpark" 
estimate  for  the  geometric  factor  of  a rectangular  collimator. 

We  used,  however,  the  exact  expression  given  in  Appendix  B in 
deriving  our  final  results. 

Monte  Carlo  Program.  In  Appendix  C is  a listing  of  the 
Monte  Carlo  computer  code  used  to  analyze  instrument  response. 
The  basic  idea  is  to  trace  electrons  through  the  instrument  to 
determine  if  they  hit  the  channeltron.  In  this  manner  the 
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effect  of  design  can  immediately  be  recognized. 

Input  data  includes  aperture  geometry,  plate  radii,  arc 
length  of  plates,  nlate  voltage,  channeltron  size,  and  the 
position  of  the  aperture  and  channeltrons  relative  to  the 
plates,  A retarding  potential  can  also  be  inputed  which 
impedes  secondary  electrons  produced  in  the  plates  from 
hitting  the  channeltrons.  Provision  has  also  been  made  to 
examine  elastic  scattering  off  the  plate  surfaces.  This 
feature  will  not  be  used. 

Let  us  now  consider  an  isotropic  flux,  J,  incident  on  the 
outside  aperture.  The  number  of  incident  particles  per  unit 
area  is  given  by 

dN  = J sinS  cos0  d0  (1) 

or  if  R is  a random  number  over  the  interval  0-1  then  the 
relation 

0 = %cos-l(l-2  R)  (2) 

generates  the  appropriate  angular  distribution  of  incident 
particles . 

As  a test  we  generated  6 x 10^  particles  through  the 
high-energy  collimator  using  equation  (2).  The  theoretically 
calculated  geometric  factor  for  this  collimator  (1.8  cm  high 
X 0.2  cm  wide  x 1.15  cm  long)  is  0.0622  cm^/sr  (Appendix  B)  . 
The  "measured"  geometrical  factor  using  particle  tracing  was 
found  to  be  0.0624  cm^/sr  which  agrees  with  the  theoretical 
value  to  within  0.32%.  This  result  confirms  that  the  computer 
code  was  properly  simulating  an  incident  isotropic  flux 
through  the  collimator. 
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IV.  RESULTS 


The  Method.  Instrument  calibration  is  determined  by 
comparing  the  predicted  (Monte  Carlo)  response  curves  with 
the  measured  ones.  First,  we  look  at  the  energy  resolution 
for  a normally  incident  beam.  Then  angular  response  is 
examined  in  directions  parallel  and  perpendicular  to  the 
plates.  After  showing  substantial  agreement  with  experimental 
data,  the  Monte  Carlo  code  is  used  to  calculate  the  final 
normalizations  for  an  incident,  isotropic  flux.  Comparisons 
are  made  with  results  from  approximate  analytic  techniques. 

Channeltron  detection  efficiency  as  a function  of  incident 
electron  energy  is  also  discussed  and  included. 

Throughout  this  section  we  use  energy  channels  8 and  9 
to  determine  the  response  properties  of  the  two  head 
as  s emb lies. 

Energy  Resolution.  In  Figure  5 the  Monte-Carlo  and 
experimental  results  are  shown  for  a normally  incident 
electron  beam.  The  beam  energy  was  systematically  scanned 
across  the  channel  to  obtain  the  energy  response  curve.  It 
is  seen  from  Figure  5 that  the  Monte  Carlo  and  experimental 
results  are  in  excellent  agreement. 

Angular  Response.  Figure  6 shows  the  angular  response 
(a  ) of  the  high  energy  plates  in  the  plane  perpendicular 
to  the  plates  (R-plane).  In  this  case  the  measured  resolution 
is  somewhat  wider  than  that  predicted  by  Monte  Carlo.  This 
effect  is  probably  due  to  scattering  off  the  plate  surface. 

Figure  7 shows  a similar  response  curve  for  the  low-energy 

plates.  Observe  that  there  are  fewer  measured  particles 

than  predicted  at  large  positive  angles.  This  effect  is 

probably  from  fringe  electric  field  resulting  from  larger 

plate  separation.  See  Figure  2.  The  triangular  shape  of  the 

theoretical  response  curve  indicates  that  the  front  collimator  . 

I 
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is  the  dominant  determinant  of  the  angular  response. 

Now  we  look  at  the  angular  response  (3  ) in  the  plane 
parallel  to  the  plate  surfaces.  All  angular  response  curves 
are  determined  at  the  central  (peak)  channel  energy.  Figure 
8 shows  the  appropriate  curves  for  the  high  energy  (narrowly 
separated)  plates.  Note  that  while  the  resolution  is  in 
excellent  agreement,  there  is  a systematic  shift  of  about 
2.5°.  Because  the  resolution  is  about  8.0°,  we  do  not  ] 

consider  this  error  as  being  serious.  Figure  9 shows  similar  j 

curves  for  the  low  energy  (widely  separated)  plates.  Here 
there  is  also  a systematic  shift  of  about  1.25°.  Note  that 
for  both  angular  scans  the  measured  resolutions  for  the 
low  energy  plates  are  less  than  the  predicted  ones.  This  is 
opposite  to  what  one  would  normally  expect  and  indicates  a 
slight  systematic  error. 


Final  Normalizations.  The  final  normalizations  were 
determined  using  the  Monte  Carlo  program.  As  described 
previously,  an  isotropic  flux  is  taken  incident  on  the  front 
aperture.  One  thousand  electrons  are  traced  through  the 
plates  at  each  energy.  Figure  10  and  11  show  the  results  for 
channels  8 and  9 respectively.  The  ordinate  represents  the 
percentage  of  electrons  that  are  detected  after  leaving  the 
collimator  (i.e.  the  transmission  efficiency  through  the  plates). 
The  integral  of  this  curve  times  the  aperture  G-factor  gives 
the  final  normalizations. 

For  an  incident  flux  J the  count  rate  (CR)  is  given  by 

JG  = CR 

where  G is  the  normalization  factor  in  cm^/sr/ev.  The  G- 
factor  can  be  approximated  by 
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G = An  A E 

A = area  of  entrance  aperture 
n = acceptance  solid-angle 
E = energy  bandwidth 
Dividing  by  the  center  energy  E we  have 

H = G/E  = An  ^ . 

E 

Now  since  A,  n and  A E/E  are  constant  for  each  set  of  plates, 
then  this  one  number  H characterizes  the  normalization  for 
all  channels.  For  the  high-energy  (narrowly  separated  plates) 
channels  H = 1.30  x 10"4  cm^/sr.  For  the  low-energy  (widely- 
separated  plates)  H = 4.3  x 10“^  cm^/sr. 

We  now  compare  this  result  with  an  approximate  analytic 
approach.  In  this  approach  it  is  assumed  that  the  energy  and 
angular  response  curves  are  independent  gaussian  distributions. 
The  final  normalization  is  the  integral  product  of  these 
distributions.  The  result  is 

H = 1.20  AG  Rq  AaAB 

where  A is  the  aperture  area,  Rq  is  the  channeltron  detection 
efficiency  (at  E 1 keV,  g = i),  Rq  is  AE/E,  Aa  and  A6 
are  the  FWHM  (full-width  half  maximum)  of  the  angular  response 
curves  . 

Using  the  measured  results  tabulated  in  Table  1,  this 

- 5 2 

approximation  gives  us  H = 6.8  x 10  cm  /sr  for  the  high 
energy  plates  and  H = 5.55  x 10"^  cm^/sr  for  the  low  energy 
plates.  Note,  however,  that  energy  and  angular  resolutions 
are  not  independent  for  this  particular  electrostatic 
analyzer,  contrary  to  the  assumptions  in  the  approximate 
approach . 
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A wide  collimator  allows  more  electrons  at  various 
energies  to  reach  the  channeltrons . A narrow  collimator 
allows  fewer  electrons  to  reach  the  detectors  over  a 
narrower  energy  bandwidth.  This  effect  is  seen  by  comparing 
Figure  5 with  Figures  10  and  11.  The  energy  resolution  of 
the  high-energy  plates  is  2.3  times  greater  for  an  isotropic 
flux  than  for  a normally  incident  flux.  For  the  low  energy 
plates  it  is  1.5  times  greater.  We  interpret  this  difference 
to  be  due  to  the  collimator  in  front  of  the  high  energy  plates 
having  a G-factor  16  times  that  of  the  one  in  front  of  the 
low-energy  plates.  It  is  for  this  reason  we  consider  the 
Monte  Carlo  normalization  as  more  closely  reflecting  realistic 
fluxes  . 

The  final  normalizatin  values  for  each  energy  channel 
as  given  in  Table  2.  Dividing  the  count  rate  by  the  approx- 
imate normalization  gives  the  equivalent  flux  in  electrons 
per  cm^/sr/ev/sec . The  channeltron  efficiencies  were 
obtained  from  Archuleter  and  Deforest  (1971).  For  1 keV£ 

E<  50  MeV,  e=  1.0  - 2. 0/(3. 0 + 6.5/(E  - 0.5)  + 30.0/ 

(E  - 0.5)2)  where  E is  expressed  in  keV.  For  10  ev£E£70  eV 
e was  taken  to  be 

E = 0 . 10  E® • (E  in  ev) 

Note  that  the  midpoint  energies  in  Figures  10  and  11  ar 
not  the  same  as  in  Table  2,  which  are  the  final  calibrated 
values.  This  does  not  affect  H,  however,  which  is  energy- 
independent.  The  stated  errors  in  Table  2 for  the  various 
energy  channels  are  a best  estimate  from  electronic  and 
experimental  uncertainties. 
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TABLE  1 


Measured  Energy  and  Angular  Resolutions  With  a 
Monoenerget i c and  Unidirectional  Electron  Beam 


Channels  1-8 


Channels  9-16 


1 . 6±0 . 2° 


3 . 7±0 . 7° 


8 . 0±0 . 8‘ 


4 . 75±0 , 08° 


AE/E(R^) 

G (Aperture) 


4 . 0 ± 0 . 4 % 

6.22x10  ^cm^-sr 


7 . 2 ± 0 . 4 % 

- 3 2 

3 . 929x10  cm  -sr 


TABLE  2 


SSJ/3 

Normal i z at  ions 

H = 1 . 

30xl0~^cm^-sr 

G = HEe 

Channel 

Center 

Channe 1 tron 

Normalization  (G) 

Number 

Energy  (eV) 

Efficiency  e 

HEe(cm  -sr-eV) 

1 

20,000±3% 

0 ,405 

1.05x10° 

2 

13 , 700±3% 

0 .43 

7 . 66x10*  ^ 

3 

8,990±3% 

0.48 

5 . 59x10*  ^ 

4 

5 . 500±3% 

0.56 

4.00x10*^ 

5 

3 . 790  + 3% 

0 . 66 

3. 25x10*^ 

6 

2 ,290  + 3% 

0 . 83 

2.47x10'^ 

7 

1,590±5% 

0,94 

1 .94x10'^ 

8 

1,060+6% 

0.99 

1 .36x10*^ 

H = 4 . 

3x10  cm  -sr 

9 

1 ,045±3% 

1 . 00 

4 .49x10*^ 

10 

661+3% 

1 . 00 

2 . 84x10*^ 

11 

434±3% 

1 . 00 

1 . 86x10*^ 

12 

264±3% 

1 .00 

1.13x10*^ 

13 

183±3% 

1 .00 

7 . 86x10*^ 

14 

1 10±3% 

1 . 00 

4 .73x10*^ 

15 

77±5% 

1 . 00 

3.31x10*^ 

16 

5 1+6% 

. 76 

1.67x10*^ 
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ENERGY  RESPONSE  - ISOTROPIC  INCIDENT  FLUX 
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1000  PARTICLES  TRACED  FOR  EACH  ENERGY 


ENERGY  RESPONSE  - ISOTROPIC  INCIDENT  FLUX 
CHANNEL  9 - MONTE  CARLO  PROGRAM 


APPENDIX  A 


ELECTROSTATIC  ANALYZER  (ESA) 
SSJ/3  OR  GFE3 

Summary  of  Characteristics 


Particles  Detected 


Electrons 


Detectors  Used 

Method  of  Energy  Analysis 
Number  of  Energy  Bins 


2 Channeltron  Electron  Multipliers  for 
each  of  the  two  sets  of  ESA  plates 

Voltage  Stepping  on  ESA  plates 

8 for  each  set  of  ESA  plates;  total  of  16 


Energy  Bins 

Large  Plates 

Sma 1 1 Plates 

20.000  keV 

1.000  keV 

13.037 

,652 

8 .498 

.425 

5 .539 

.277 

3.611 

.181 

2 . 354 

.118 

1 .534 

.077 

1 . 000 

.050 

Acceptance  Angles 

Large  ESA: 

1.6° 

8 . 0° 

FWHM 

FWHM 

across  the 
along  the 

apertures 

apertures 

Sma 1 1 ESA : 

3 . 7° 

4 . 8° 

FWHM 

FWHM 

across  the 
along  the 

apertures 

apertures 

Normalization  Constants 

Large  ESA: 

1 . SOxlO'"^ 

2 ^ 

cm  -ster. 

Small  ESA: 

- 3 2 

4.3x10  cm  -ster . 

Large  ESA: 

4.0% 

E 

Small  ESA : 

7 . 2% 

Data  Rate  One  complete  Spectrum  per  second 

Dwell  Time  at  Each  Energy  98  ms. 

Level 


Digital  Data  Format 


(9  bits  per  channel)  x (16  channels)  = 

144  bits 


Note:  Each  9 bits  in  every  channel  are  in  logarithmic  form,  the  five 

least  significant  bits  being  the  mantissa  and  the  remaining  four 
the  exponent.  This  number  is  converted  to  decimal  form 
according  to  the  following  relationship: 

N = 2^(x+32)  - 33 


APPENDIX  B 


The  G-factor  for  a rectangular  shaped  collimator  can 
be  calculated  exactly  (Rothwell  and  Moomey,  1972). 

With  the  front  and  rear  apertures  having  the  dimensions 
(^2’  ^2^  respect  ive  ly  ^ one  can  define 

a = a^,  a^  width 

6 = h(.b^+h^)  b^,  height 

Y = J^Ca^-a^) 

6 = ^(bj-b^) 

If  the  collimator  has  a length  L then 


G 


r,  2 2 .2  , 2 2 .2  -1 

L +a  +6  L +Y  +B 

In  L 2 2 .2  * 2 2 ^,2 

'-L  +a  +B  L +Y  +0  -J 

+ F(a.e)  + F(y,6)  - F(a,6)  - F(y,6) 


where  F (a , S) 

F(a,e) 


F(a,6)  + F(B,a) 
26(L^  + a^)"^  tan 


r 3 


L(L^  + a ) " J 


This  expression  is  consistent  with  that  given  by  Willis  and 
Thomas  (1972)  and  by  Sullivan  (1971)  (with  a corrected 
typographical  error). 
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APPENDIX  C 

Listing  of  computer  code  that  was  used  for  the  Monte-Carlo  computations. 


PROGRAM  tSAPLTSdNFUT, OUTPUT) 

OlMCNSrON  OAT (24)  |E FF ( 200) f EE (2  00 ) , EG (20 0 ) , AMAT  (10,10),S 

♦ UMKIO) 

DIMENSION  THH(IOO), AX(IOO) ,THM( 100) 

COMMON/Fl/EFF,EG ,nELE ,N5,GINT 
R2 (X) =A*X4B 

P( X,Y, AL)=2tP* V»SQfiT ( AL**2+X**2) *ATAN(YXSQRT( AL**2+X**2) ) 
N5=60 
N7=N5/6 
Nq=5 

Wll=3fi 
N12=Nll/2 
AMF=l,€6E-27 
Q0VMP=9.638E  07 
AME=q.lE-31 
Q0ME=1.7582E  11 
AA=  180./  2.141596 
CNNTfl=0.0 
ELEC=8ME  lECTRON 

1 continue 

CALL  RANSET(TIME(CUM)) 

C)'NTA  = 0. 0 

C IF  11=1  nc  BACKGROUND  CALCULATION  EASE  ON  "LINE  OF  SIGHT- 
C TJ.NE.  0,  automatically  ASSUMES  PARAMETERS  FOR  L.  FRANKS  E 

♦ SA 

C IJ=2,  SIMULATES  ELECTRON  GUN  SCAN  OF  FRANKS  ESA  AS  SHOWN  IN  F 
♦IG.  2 OF 
C HIS  PAPER 

t C 11=1,  DOFS  LINE  OF  SIGHT  BACKGROUND  CALCULATION 

C K1=0,  NO  ELASTIC  SCATTERING  OFF  CF  ESA  PLATES 

C Kl=l,  MORE  THAN  1 SCATTERING  OFF  OF  ESA  PLATES  PEP  ELECTRON  A 
♦S  DEFINED 

C BY  K2  (INNER  PLATE) 

C Kl=2,  1/K2X100  PER  CENT  SCATTERING  OFF  OF  ESA  PLATES  (INNEF  P 
♦LATES) 

C Kl=3,  SCATTERING  OFF  OF  OUTER  ESA  PLATE  AS  DEFINED  BY  K9 
C Kl=4, SCATTERING  OFF  OF  OUTER  tSA  PLATE  AS  DEFINED  BY  1/K9X100 
C Kl=5, SCATTERING  OFF  OF  BOTH  PLATES, K2  TIMES  INNER  PLATE, K9  T 
♦IHES  CFF 
C ATE 

C K1=6,CATTERING  off  OF  BOTH  PLATES1/K2  INNER, 1/K9  TIMES  OUTER 
♦PLATE 

C IF  IK=0  INCIDENT  FLUX  IS  ISOTROPIC 
C IK  = l INCIDENT  FLUX  AT  INCIDENT  ANGLE  ANGE 
C IF  IK=2  SCAN  IN  R-PLANE  IN  STEPS  OF  ANGE  DEGREES 
C IF  IK=3  SCAN  IN  Z-PLANE  IN  STEPS  OF  ANGE  DEGREES 
C NOTE—  IF  IK  GE  2 INPUT  DATA  PACK  CHANGES 

READ  4,II,IJ,K1,K2,K9,ANGA,ANG2, ANGD,IK,ANGE 
4 FOPMAT(5I5,3F7.2, I5,F7.2) 

PFINT  1900, IK, ANGE 

! 19CC  FORMAT  (IX, IS, F7. 2) 


10 

11 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 
2 6 
J7 
28 

29 

30 

31 

32 

32 

33 

33 

34 

35 

36 

37 
37 
3ft 
39 
3a 


4 2 

42 

43 

44 

44 

45 
4 6 
47 

46 

49 

50 

51 

52 

53 


IF<K1.EO.O)  K7  = l 
IF(K1,£0,1)  K7=l 
IF(K1.EQ.2»  K7=-l 
IF(K1.EQ.3I  K7=2 
IFCKl.EO.M  K7=3 
IF(K1,E0.5»  K7=2 
IF(K1.ECI.6>  K7  = 3 
IF(K1.EQ.0I  K8=l 
IF(Kl.EQ.l)  K8=2 
IF(K1.EQ.2)  K8=3 
IF<K1.EQ.3)  K8=l 
IF<K1.E0.«»)  K8  = l 
IF(K1,EQ.5»  K8=2 
IF(Kl.eQ.6)  K8=3 
IF(K1.EQ.7)  K7=2 

IF<Kl.EO,7l  K8=3 
1FCK1.EQ.8)  K7  = 3 
IF(K1,E0.8>  K8=2 
READ  3,  COATd)  ,1=1,24) 

3 FORMAT! <12A6) > 

REAn2, FART, POST 

2 FOPHAT(A8,F3.0) 

PRINT  5,  WATd)  ,1  = 1,  24) 

F FORmaTCIH,  12A6/12fi6l 

IFdI.EQ.l)  PRINT  10 

10  fopnat(ix,*9ack3round  calculation*) 

IFdK.EQ.l)  PRINT  12,ANGE 

12  FORKATIIX,*  CONST.  ANGLE  SORT ,A NGLE=* , FT . 2 ,♦  CEG.*) 

IF  (oarT.EQ.E  IEC»  AMP  = AME 
IF<nART,EQ.ELEC)  QOVNF  =-OOME 
PRINT  6, PART 

f.  FCRRAT!  1X,*TRACING*,3X,A8/) 

IFCDART.EQ.ELEC.ANC.IJ.NE.  G)  POST  = 160. 

IF(IJ.NE.O)  PRINT  16 

15  FOFNATdX,*  FRANKS  ESA  ♦) 

IFdJ.EQ.2)  PRINT  13 

18  FORMAT  ( IX, •ELECTRON  GUN  FLUX  = 2.0E  06  ELEC/CMSG/SEC*) 

IF(K1,NE.0)  PRINT  17,K1,K2,K9 

17  FORMATdX, •ELASTIC  SCATTERING  OFF  OF  ESA  PLATES  ,*, 

♦•<!=•, 15, 

A • K2=*,I5/1X,*IF  Kl=l  K2  SCATTERINGS  PER  FART.,  IFK1  = 2 

♦ EVERY  OT 

eHEP  K2  PART,  MAY  SCATTER* ,*K ?=♦, 15) 

IFCKl.NE.O)  PRINT  21,ANG2 

IF(K1,NE.0)  PRINT  19,ANGA 

IF(K1.N£.0)  PRINT  23,ANGD 

19  FOFMATdX,*  MAX  SCATTECING  ANGLE  OFF  OF  OUTER  pLATE=*,F7. 

*2,*  DEG*) 

21  FOFMATdX,*  MAX  SCATTERING  ANGLE  OFF  OF  INNER  PLArE  = *,F7. 

♦2,*  DEG*) 

23  FC°MATdX,*MAX.  ANGLE  FOP  INELASTIC  SCATTERING  CPF  OF  OU 

♦TEF  PLATE 


'I 
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best  AVAIUEIE  COPY 


A»*,F7.2)  101 

C SFT-UF  APERTURE  AND  CALCULATE  ITS  G- FACTO*?  l'?2 

»EAO  e,Al,Pl,A2^B2,ALl,AL2,  AL  3 , \/OLT  1 , K5  103 

6 F0P.f'AT(7F7.3 ,2F7.3)  lOL 

PRINT  ll,Al,Bl,A2,e2 ,AL1,AL2,AL3 ,V0LT1 ,R5,P0ST  lOE 

11  FOR“AT{lX,*APERTURfcDIhENS10NS*/lX,*FP0KT  WIOTH=*  ,F7.  3 ,*  C ICC 

♦ M*,3^(,*FR  lOe 

10NTHT=*,F7. 3f*  CM*,3X,  *BACK  W IOTH=» , F 7. 3 , * CN*,3X,*BACK  107 

♦ HT=*,F7.3  1C  7 

2,*CM*/1X  ,*LENGTH  of  A PEP. . =•  , F7 . 3 ,♦  CM*,3X,  *OIST  BTMN  APE  108 
♦R  AND  FLA  108 

3TFS=»,F7.3,*  CM*/1X, *OIST  BTWN  PLATES  ANO  CHANN , =* ,F 7. 3, ♦ 109 

♦ CF*,3X,*  1P9 

4B1AS  VOLTAGE  IN  FPCNT  OF  DET ECT, =* , F7. 1 , * VOLTS* /IX, •CF AN  110 
♦NELTRCF  R 110 

1A0IUS=*,F7.2,*  CM*,*  POST  ACCELERATI ON= * fFS  . 0 , *VOL T S*/ ) 111 

ALFHA=(A14-A2I/2.  0 112 

BETA=  {ei+B2» /2.0  113 

GAH=<A1-A2I /2.  0 llu 

OELTA  = (01-02)/:2.O  115 

G1  = AL1**2*AL0G(  ( AL  1*  * 2+ALPHA**  2*  DELTA*  *2  ) * ( AL  1**  2+&AM**  2*  HE 

*BETA**2)/  llfi 

K AL1**2*ALPHA**2*B£TA**2>/<AL1** 2 ♦GAM**2*0EL TA**2)  I 117 

G2=F(  ALPHA,BET  A,AL1)  fF  (BETA,  ALPHA,  ALl)  t-F  ( G AM,  CELT  A,  ALl  > *F  113 

♦ (DELTA, GA  118 

1M,AL1)-F  (ALPHA  , DELTA  , ALl ) -F(  OELT  A,  ALPH A , AL  1 ) -F  (G AM,  BE TA  , A 119 
♦L1I-F(PET  119 

2A, GAM, ALl)  120 

G=G1+G2  121 

PRINT  13,5  i;?2 

13  FORHATC  1X,*G-FACTCR  FOR  APERTUR£=*  ,E1  C.4  ,♦  CMSO-SR*/)  123 

PRINT  15  124 

IF  F0FMAT(1X,*ESA  PLATE  DIMENSIONS*)  125 

C R2A  = RACIUS  AT  APERTURE  12f. 

C R2P=RACIUS  AT  EXIT  127 

PEAO  7,R1,R2A,R28,THE,V0LT,H,T0T,R9  129 

7 F0FMAT(4F7.3,  E10.2,  F7 . 3 ,E 10  .4 , F7 . 2)  129 

PFINT  9,Rl,R2A,R2e,TME,V0LT,H,T0T,R9  130 

9 FOf'MAT(lX,*Z£RO  POTENTIAL  RA  01  LS=*,  F7.  3,  * CM*  , 3X , *RADI  US  A 131 

♦T  FOTENTI  131 

lAL  V AT  AP£RTURE=*  ,F7.3,*  CM*  , 3X  ,*RAOI US  AT  POTENTIAL  V A 132 

♦T  EXIT=*,  132 

2F7.3,*  CM*/3X,*MAX.  ANGLE  OF  PLA TE S=*, F7 . 2 ,♦  DEGS*,3X,*V0  133 

♦ LTAGE  ON  133 

3 PLATeS=*,lPE10. 2,*  VOLTS*  /3X,*PLATE  HEI GHT=* , 0 PF7. 2 , * 134 

♦CM*,  * N 134 

LUMEEP  CF  PARTICLES  TRACED  FOR  EACH  ENE FG Y=* ♦ IFE 1 0. 2/ 1 X , *0  135 

♦ISTANCE  0 135 

5F  CHANNELTRON  FROM  CENTER  OF  ESA  AXIS=  * , IPEl 0 . 2 , * CM*/)  1.16 

IF(PART,EQ,E  lEC)  VCLT--VOLT  137 

THE=THE/AA  136 


(»«“>• 
1 I 


OE-02 

& NGB“ AKoA/ A A 
aNC.G  = ANG2/AA 
ANGr=ANG^/AA 
AhG‘^=ANGe/ AA 
Rl=oi*l,0E-02 
R2A=R2A*1.0€-02 
R2B=P2'?*1.0E-02 
P5  = i’5*1.0E-02 
R9=R9*1.3E-02 
9=R2A 

A=  CR23-R2A) /THE 
WIG1= AES(R2A-P1I 
RR=(R1*R2B) / 2, 0 
A1  = A1*1.  CE-02 

31  = 91^1.  (E-0  2 
A2  = A2*l.GE-0  2 

32  = 82*1.  jE-G2 
AL1=AL1*1. OE-02 
AL2=AL2*1.0E-02 
AL3= AL3*1. OE-02 
cr,=  p,D-Ai/2, 0 
R7=rp-A2/2. 0 
R8=RP  + A2/2,  0 
IFdJ.KE.O)  Rb  = 1.2‘*E-01 
IFdJ.NE.O)  R7=R& 

IF(IJ.NE.O)  S8=R6-*-A2 
IF{TJ.E0.2)  THY  = -!5.0/AA 
OELV=l.a/AA 

thw=thy-cely 

TH>=THE/2.0 
P3  = 02<THX) 

GO  TO  14 

24  IF(IK.GE.2)  READ  26, VOLT 
IF(IK.GE.2)  PRINT  29, VOLT 
26  FCPHAT(  IPEIO.  2) 

29  FORPATCIX,*  VOLTAGE  ON  PL  AT ES=* , IPEIO . 2 , * VCLTS*» 
IFrPAPT.EQ.ELEC)  VClT  =-V0LT 
CONTIMJE 

C USING  HKS  SYSTEH  — ENERGY  IK  JOULES 

EE(1)=VOLT/2.0/ALOG(R3/R1)*1,6E-19 
EE (1)=ABS(EE(1)) 

CONST=GOVMP* VOLT 
EA=EE(  1)/1.&E-19 
PRINT  31, EA 

31  FOP^AT  ( 1X,*HE  AN  ENERGY=*  , IPEl  C.  2,  *E  V*  ) 

IF(IK.LT.2»  GO  TO  34 
THH( 1)=-ANG£*N12 
DO  33  1=1, Nil 
THH(I+1)  = THH(I)  -fAKGE 
THM(I»=THH{I) *AA 

BtSl  AVAllABlt  COPY 


139 
14  0 
H*1 

142 

143 

144 

145 


146 

147 
149 
149 
IfJO 

151 
1F2 

153 

154 

155 

156 
15  7 
159 

159 

160 
161 
162 
153 

164 

165 

166 

157 
168 

169 

170 

171 

172 

173 

174 

175 

176 
IY7 

178 

179 
18  0 
13  1 

152 
19  3 
194 
18  5 

18  t 

18  t 
18  9 
169 


, : 


34 


Ax<i)=o.o  ig- 

33  CCNTINUE  ISl 

THm(N11«-1I=TMH(N11«-1)*AA  19? 

GO  TC  36  193 

7k  COHTINUE  194 

OELE=0.005*eECl)  195 

EF(1) =0.9*EE(1I  196 

DO  35  K = 1,N5  197 

EE<K*ll  = EE(K)^OELE  19*^ 

35  CONTINUE  199 

36  CONTINUE  200 

DO  22  J=l,10  201 

no  22  K=l,10  20? 

AMAT<K,J»=0. 0 203 

22  continue  204 

25  CONTINUE  205 

THW=THK*CELY  206 

THW=-4,0/AA  207 

00  600  K=1,N5  208 

K13=l  2C9 

37  CONTINUE  210 

K4=0  211 

N4=0  212 

N6=0  213 

K5=n  214 

CNT=0.C  215 

COUNT=0.0  216 

CNTA=0.0  217 

C OELT=TIrtE  INTERVAL  FOR  EACH  STEP  218 

GO  TC  40  219 

38  IF(IK,€E.2)  GO  TC  202  220 

40  CONTINUE  221 

K3=n  222 

N3=0  223 

K11=0  224 

K12=0  225 

EX=EEIK)  226 

IF(IK.GE.2)  EECK)=EA*1.6E-19  227 

V=SQRT(2,»E£{K» /ANF)  226 

C PULL  OUT  DEFOREST  CARDS  229 

P=R6  + A1-*RANF<0UN)  230 

Z1=31*PANFC0UM)  231 

C CHOOSE  T1  SUCH  THAT  FOR  JCFLUX)  = CONST.  T1  IS  PICKED  SIN(T1»  232 
♦COS(Tl)  232 

0 TIRES,,  SIN  FOR  SOLID  ANGLE  COS  FOR  PROJECTED  AREA  233 

Ti=ACOS(1.0-2.0*RANF(DUM»)/2.  C 23- 

T2=6. 2831*RANF (OUM)  235 

IFdK.EQ.ll  T1=AN6E  236 

IF<IK.EQ.l)  T2=0.0  237 

IFdK.EO.2)  T1  = THH(K13)  238 

IFCIK.EQ.2I  T2=0.0  230 
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IF(IK.EQ.3»  T1=THH(K13>  2%0 

IF(IK.E.1.3I  T2=l.f707  241 

IF(IJ,£0.2I  ri-=A0S(THM)  242 

IFII  J,£0.2.  ANO.THW.GE.  C.  C»  12=180. /Aft  243 

IF(IJ.EQ.2.ftNO,TrtW,LT  .0.  0)  T2=0.  0 244 

Y1=0  24y 

V»<=V*SIN  (Tl)  ♦OOS<T2>  246 

^TT=V*COS<Tll  247 

VZ=V»SIN(T1)*SIM<T2)  248 

DT=AL1/VTT  249 

r=R4.VR*OT  2»;o 

Zl=Zl  + (/Z*DT  251 

IF(Z1. 6T.B2.0R.Z1.lt. 0.0)  GO  TO  38  252 

IF(R.GT.R8.0R.R.LT.R7)  60  TO  38  253 

CNTft=CNTA-H.O  254 

DT=AL2/VTT  255 

R=F+VR*DT  256 

Z1=Z1  + VZ*DT  257 

IF(R.GE.R2ft.0R.R.LE.Rl)  GO  TO  202  258 

ANGf'  = VTT*R  259 

TT1=0.C  260 

R3=R2<TT1)  231 

CONTT=CONST/ftLOG  (R3/R1)  252 

CONSS  = CO»^ST  263 

IFdl.EO.ll  CONSS  = C.G  264 

DELTT=C0NTT/{R1+R3)*2  265 

DrLTT=3.  IE-01*V/OELTT  266 

DELT=3.1416^THE'^(Rl«-R3)/2./\//2.0E  02  267 

IF(DELTT.LE.OELT)  CELT=0ELTT  263 

GO  TO  20  269 

27  continue  27  0 

IF|<12.GT.O.OR.K11.GT.OI  R=R*s>R*DELT  271 

GOTO  28  272 

20  CONTINUE  273 

R3=P2<TT1)  274 

CONS2=  {VTT»*2-CONSS/ ALCG<P.3/R1))  /R  275 

r=p+VR*C£LT+CONS2*OELT*»2/2. 0 276 

VP  = VR»-  C0NS2*DELT  277 

28  CONTINUE  278 

VTT=aNG''/R  279 

TTl=TTl+rfTT*OELT/R  280 

Z1=Z1*  OELT^VZ  261 

IF  (Zl.LE.O.O.OR.Zl.GE.H)  GOTO  2G2  282 

IF<TT1  .GE.THE)  GO  TO  400  293 

IF  (R.GT.fil.ftND.  R.LT.P3I  GO  TO  1 700  28^ 

GO  TO  1800  295 

1700  CONTINUE  296 

<ll=r  267 

K12=0  283 

GO  TO  20  289 

180  0 continue  29'' 
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1 ,1 1 . 1 IIIIPL^^WIP 


C t<l  = l,2 SCATTER  OFF  OF  INNER  SURFACE  ONLY  291 

C *^1=3, SCATTER  OFF  OF  OUTER  SURFACE  ONLY  292 

IF(R.GE.RJ)  GO  TO  1300  293 

IF(R.LE.Rl)  30  TO  1400  294 

GO  TO  202  295 

130C  CONTIKUE  296 

IF(Kll.GT.O)  GO  TO  27  297 

K11=K11+1  29ft 

GO  TO  (2C2,  1100,  1200)  ,t<7  299 

ICOC  CCMTIKL'E  30  0 

IF  (K12.GT,  0)  GOTO  27  30  1 

K12=t<12+1  302 

GO  TO  (202,900 , 10  00)  , KB  303 

900  continue  304 

N3=N341  305 

IF(N3.GT.K9)  GO  TO  2C 2 306 

VO=AeS(VR)  307 

THF  = ATAN(VD/VTT)  305 

IF (APS  (THP) .GT .ANGC)  GO  TO  1500  309 

VF--VR  310 

GO  TO  27  311 

1000  CONTINUE  312 

VO=A-3S(VR)  313 

7HF=ATAN(VD/VTT)  314 

Nt=N6+l  315 

rF(A5S(THP> .GT. ANGC)  GO  TO  1500  31F 

VF=-VR  317 

GO  TO  2T1  3ie 

HOC  continue  319 

VO=A9S(VR)  320 

K3=K341  321 

IF(K3.GT.K2)  GO  TO  2C 2 322 

THP=  ATAN(i/R/yTT)  323 

IF(ARS(THPI .GT.AN6C)  60  To  202  324 

IF  (APS  (THP)  ,GT  .ANG6)  GC  TO  1600  325 

VF=-VR  32b 

GO  TO  27  327 

1200  CONTINUE  328 

V0=A9S(VR)  329 

K4=K4+1  330 

THP=  ATAN(i/R/yTT)  331 

IF(A9S(THP) .GT. ANGC)  GO  TO  202  332 

IF(A0S(THP) .GT.ANGfi)  GO  TO  1600  333 

VR=-VR  334 

GO  TO  200  335 

1500  CONTINUE  336 

EX  = EX  *(1,0-RANF(DUH,))  **(1.0/1.35)  337 

VA=SQPT (2.C*EX/AMP)  33£ 

VR=VA*SIN(THP)  339 

VTT=VA*CCS(THP)  340 

IF(K8.EQ.3)  GO  TO  201  341 


T>^f- 
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GO  TO  27 

342 

160  C 

CONTINUE 

343 

EX  = EX  * (1.0-RANF  (OUM) )•* (1.0/1.35) 

344 

VA=SQRT(2.0*EX/AMP) 

345 

VR=- WA*SINfTHP) 

346 

VTT  = VA*COS(  THO) 

347 

IF(K7.EQ.3)  GO  TO  200 

348 

GO  TO  27 

349 

200 

CONTINUE 

350 

IAe=MC0(K4,  <21 

351 

IF(IA8.£Q.O)  <4=0 

352 

<5=<4 

353 

IFdAP.EQ.O)  GO  TO  27 

354 

GO  TO  202 

355 

201 

CONTINUE 

356 

IAC=NCC(N6,<9) 

357 

IFdAC.EQ.O)  N6=0 

358 

Nt.  = N6 

359 

IFdAC.EQ.O)  GO  TO  27 

3&' 

2C2 

CONTINUE 

3S1 

CNT=CKT-H 

362 

IF (CNT.GE.TOT)  GO  TO  500 

353 

GO  TO  40 

384 

4CC 

CONTINUE 

365 

C SET 

UP  matrix  AT  EXIT  APERTURE  OF  PLATES 

3 36 

00  t20  N=1,N9 

367 

IF(Z1  .G£.(N -1)*H/ N9.AN0.Z1.lt. N*H/N9)  N1  = N 

368 

<!2r 

CONTINUE 

359 

00  430  N=1,N9 

37  G 

RRR1  = R1+  (N-1)  ♦WI01/N9 

371 

RRF2=R1*N*WIC1/N9 

372 

IF<R.GE.RRR1.  AN0.R.LT.KRR2)  N2=N 

373 

43C 

CCN'»’INUE 

37*. 

AMAT(N1,  N2)  =AMAT(  N1,N2)  «■!.  0 

375 

c oc 

WE  HIT  THE  CHANNELTRON? 

37  6 

RR=(Pl-*-R2B)  /2.  0 

577 

C FULL  OEFOFEST  ESA  CARDS 

378 

0IS2= AL3 

379 

0T=0IS2/ NTT 

3 8'* 

XX=0. 0 

381 

</XX=  VTT 

382 

VT=V0LT1 

383 

OTT=0 . l»OT 

384 

TT=C.0 

385 

ACC=0CVMP*VT/0IS2 

366 

1^60 

CONTINUE 

387 

TT=Tf4CTT 

338 

XX=XX  + VXX*0TT-ACC»0TT**2/2.  C 

389 

VXX=</XX-ACC*OTT 

39  0 

C DC 

ES  LOW-ENERGY  ELECTRON  STOP 

391 

IF(VXX.LE.O.O)  GO  TO  202 

392 

38 


BKf  COPY 


IF<XX,LT.0IS2)  50  10  1*60  393 

«=H-VP*TT  39*4 

Z1=Z1+V2*TT  395 

IF(IJ.NF.O)  P.5  = 1.77E-  03  396 

IF<IJ.KE.O)  R9=-RF  397 

R4=a9S (RR-R9-R)  398 

C PULL  OUT  CEFOREST  C4RDS  399 

IF(P.LT.(RR-A2/2. 0)  .0R.R.5T  .{RR  + A2/?.  0)  ) PO  TO  202  ««  0 0 

IF(71.0T.02.0R.Z1.lT.  C.  0)  GO  TO  202  901 

P«)A  = SGRT  (R4*»2*-(  Zl- 92/2 . 0-R5)  ♦♦  2 ) t0  2 

Pi*9=SQRT(P.9*'^  2<-(  Zl-9  2/2  .C+R5)  ♦♦2)  9'''  3 

IF(P«tA,5T.R5.AND.RV9.GT  ,P.F)  GO  TO  202  <!  0 L 

IF(TK.GE.2»  AX(K13)  =AX(  K13)  -t-l  .0  *.0  5 

nOUNT  = COUNT«-l  **0  6 

GO  TO  202  i;"7 

5:r  CONTINUE  0 c 

EFF(t<)  =COUNT/CMT  i*09 

IF(IK,G=’.2)  AX(K13)  =AX(K13> /CNT  410 

EG<K)  =EE(K) /I.  6E-19  411 

CNMT A=CNNTA+CNTA  412 

PRINT  5 C7,EG<  K)  , CNT  A ,CNNTA  413 

507  FO'^MATdX,*  AT  ENERGY  = ♦,1PE10.2,»  EV  NUf-^eEP  OF  PA  1 1 4 

♦RTICLES  414 

1 THPQUGH  APERTURE=  •,0PF10.2,  * RUNNING  TOTAL-  »,0PF10.2)  Rl*^ 
IF(IT.EO.l)  GO  TO  630  4i€ 

£99  = EG(K) /I.  OE  0 3 41  7 

EFI=1.0  418 

TENP=EG(K)  419 

EG  (K)=EG(K)-/T+POST  *120 

IF(PART. EO. ELEC .AND .EG(K) .GE. 1. CE  03. AnD.£G( K) . LE.5. OL  0 421 

4-6J  421 

1 FFI  = 1.0-(  2 .0/(  3.  0fF..F/ (E99  -0  .5) +30  ./ {£99  -0.5)**3))  422 

IF  (PART.  EQ.ELEC.  AND.  £G(K)  .GE.l  C.O.AND.  EG(K)  .LE.7  0.)  £FI=  0 423 

+.10*EG(K)  423 

1**0.515  424 

IF  (PA  FT.  EQ.ELEC.  AND.  EG(<  ) .GE  .2  00  . AND.  I J.  NE . 0)  EF I = 8. . 1 9 9/ EG  4 25 

+(K)+’0.4l  425 

A7  42fi 

IF  (PART.EO.E  LEC.ANC.  EG(K)  .LT.  2 CO  .AND.IJ.NE.  0)  EFI=5 .9E-04*'  42  7 

+ EG<'<)**1.  427 

A3e  428 

EG<K)=TEHP  429 

IF(IK.GE.2)  AX{K13)  =AX(K13) ’EFI  430 

IF(IK.GE.2)  GO  TO  440  431 

EFF<K>*=EFF(K)*EFI  432 

IF(IJ,EQ.2)  EFF(K» =E FF{K) *2. OE  06  433 

IF(  PCD  (K,N7)  .NE . C ) GO  TO  440  ‘•34 

E1=EF(K-N7+1»/1.6E-19  435 

E2  = EE  (K)/l.ftE-19  43f. 

AU7  = t’7  437 

00  435  N=1,N9  43^ 
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SU»«1(N)=0.0  439 

no  435  1=1, N9  440 

SlimCNJ=SUHl(M)fAM«T(  I,N»  441 

435  CONTINUE  442 

flNe=N7*TOT  443 

PRINT  43  3,E1,E2,  ((flMaT(Nl,N2»  ,N2  = 1,N9I  ,N1  = 1,N9»  , (SU^l  (t.)  , 4ut. 
♦N=1,N9),  444 

1AM8  445 

433  FORMAT  (IX, •EFFICIENCY  MATRIX  FOR  EXIT  APERTURE  OF  ESA  FLA  446 

♦TES  FRCM  44F 

lENEPGV=*,lpriO .2,*EV» ,*TO» ,1PE 10 .2,*EV^/0PPF10 .2/5F10 ,?/5  447 

♦F10.2/5F1  447 

Z0.2/nF10.2/lX,*TOTALS*/5F10. 2, *P ARTICLES  TRACED: *, IPE 1C. 2 44 € 

• //)  44  0 

DO  438  1=1,  10  44E 

DO  438  N=1,1U  45C 

AMAT  (I,n)  = (}.Q  451 

438  CONTINUE  45  2 

440  CONTINUE  453 

IF(IK.GE.2)  K13=K1341  u54 

IF(IK.GE.2.ANO.K13. LT.Nll)  GO  TO  37  455 

1F(  IK.GE.2)  GO  TO  5 10  ‘‘SE 

GO  TO  600  457 

600  continue  458 

CALL  SIMP  459 

SIC  continue  460 

IF(IX.EC.3)  PPINT  613  461 

613  FCRMATdX, ♦ANGULAR  SCAN  IN  Z-FLANE»I  452 

IF(IK,E0.2)  PRINT  609  463 

6C9  FCPMAT  (IX,  • ANGULAR  SCAN  IN  R-FLANE*)  464 

IFCIK,6E.2»  PRINT  611,  ( THM ( K> , AX C K) , K= 1 , N 11)  465 

611  FOPMAK  (1P12E 10.  2>)  466 

IF(IK.GE.2)  GO  TO  24  467 

63r  CONTINUE  466 

thf;=thk*aa  469 

IF(IJ.E0.2)  PRINT  74  0, THP  47  0 

740  FORMAT (1 X, •ELECTRONS  FROM  GUN  INCIDENT  ON  APERTURE  AT  AN  471 

♦angle  CF  471 

740  A •,P7.2,^  OEG^)  47  2 

PRINT  51  i.73 

51  FOPMATdX, •EFFICIENCY  OF  CHANNELTRON  TO  SEE  PARTICLES^)  474 

PRINT  55,(EG(K) ,EFF(K) ,K=1,N5)  475 

55  FOPMATdX, ♦ENERGy-EFFICIENCV*/(lP10E12. 2)  ) 4^6 

G1=G»GINT  477 

PRINT  63,GINT,G1  478 

o3  FORMAT  {IX,^ INTEGRAL  OF  CHANN.  EFFIC. =• ,1PE10. 4,  • EV^,3X,*  479 

♦NORMALIZA  479 

ITION  FACTOR  FOR  ENTIRE  £S A=» , IPE 10. 4, • CMSC-SR-EV^)  480 

HH=G1/EA  481 

PRINT  67, HH  452 

67  FOPMATdX, •DEFOREST  H-FACTOR=^  ,1P£10  .4  ,•  CMSO-SR*)  483 

GO  TO  1 484 

END  465 
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SUBROUTINE  SIMP  466 

■DIMENSION  EFF(  200)  ,EG(200I  tS? 

COHNON/Fl/EFF,-cG,OEL  E,K5,GINT  48  8 

C SIMPSON  RULE  h89 

DEL=DELE/1.6E-19  491 

NftsNR-S  ^91 

N7=n5-2  ^92 

SUM1=0.0  493 

00  10  J=3,N6,2  494 

SUM1=SUM14-EFF<  J)  495 

10  CONTI^UE  49e 

SUM2=0«0  497 

DO  20  J=2,N7,2  498 

SUM2=SUM2+EFF (J)  499 

2D  CONTINUE  50  0 

GINT={EFF(1)  ♦EFF  (N5-1)  ♦•4.0»SUM24-  2.  0*SUM1)  *DEL /3 . 0 !»0l 

RETURN  5C2 

END  503 
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